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ABSTRACT 

We use statistical results from a large sample of about 500 blazars, based on two surveys, the Deep 
X-ray Radio Blazar Survey (DXRBS), nearly complete, and the RASS-Green Bank survey (RGB), to 
provide new constraints on the spectral energy distribution of blazars, particularly flat-spectrum radio 
quasars (FSRQ). This reassessment is prompted by the discovery of a population of FSRQ with spectral 
energy distribution similar to that of high-energy peaked BL Lacs. The fraction of these sources is sample 
dependent, being ~ 10% in DXRBS and ~ 30% in RGB (and reaching ~ 80% for the Einstein Medium 
Sensitivity Survey). We show that these "X-ray strong" radio quasars, which had gone undetected or 
unnoticed in previous surveys, indeed are the strong-lined counterparts of high-energy peaked BL Lacs 
and have synchrotron peak frequencies, f pC ak, much higher than "classical" FSRQ, typically in the UV 
band for DXRBS. Some of these objects may be 100 GeV - TeV emitters, as are several known BL 
Lacs with similar broadband spectra. Our large, deep, and homogeneous DXRBS sample does not show 
anti-correlations between ^ poa k and radio, broad line region, or jet power, as expected in the so-called 
"blazar sequence" scenario. However, the fact that FSRQ do not reach X-ray-to-radio flux ratios and 
^peak values as extreme as BL Lacs and the elusiveness of high i/ pea k-high-power blazars suggest that 
there might be an intrinsic, physical limit to the synchrotron peak frequency that can be reached by 
strong-lined, powerful blazars. Our findings have important implications for the study of jet formation 
and physics and its relationship to other properties of active galactic nuclei. 

Subject headings: galaxies: active — BL Lacertae objects: general — quasars: general — radiation 
mechanisms: non-thermal — radio continuum: galaxies — X-rays: galaxies 



1. A PHENOMENOLOGICAL DESCRIPTION OF THE 
BLAZAR CLASS 

Blazars are one of the most extreme classes of ac- 
tive galactic nuclei (AGN), distinguished by high lumi- 
nosity, rapid variability, high polarization, radio core- 
dominance, and apparent superluminal speeds (Urry & 
Padovani 1995). Their broad-band emission extends from 
the radio up to the gamma-ray band, and is dominated 
by non-thermal processes. The blazar class includes BL 
Lacertae objects, characterized by an almost complete 
lack of emission lines, and the flat-spectrum radio quasars 
(FSRQ), which by definition display broad, strong emis- 
sion lines. 

More than 95% of all known blazars have been discov- 
ered either in radio or X-ray surveys (Padovani & Giommi 
1995a). Previous work has shown that X-ray and ra- 
dio selection methods yield objects with somewhat differ- 
ent properties at least for BL Lacs. The energy output 
of most radio selected BL Lacs peaks in the IR/optical 
(Giommi & Padovani 1994; Padovani & Giommi 1995b, 
1996). These objects are more highly polarized (Jannuzi, 
Smith, & Elston 1994), and are more core-dominated in 
the radio (Perlman & Stocke 1993; Laurent-Muehleisen et 



al. 1993; Kollgaard et al. 1996; Rector, Stocke, & Perlman 
1999) and are now referred to as LBL (low-energy peaked 
BL Lacs). By contrast, the energy output of most X-ray 
selected BL Lacs peaks in the UV/X-ray; these objects, 
which are referred to as HBL (high-energy peaked BL 
Lacs) , are less polarized and generally have preferred posi- 
tion angles of polarization in the optical (Jannuzi, Smith, 
& Elston 1994), and are less core-dominant in the radio 
(Perlman & Stocke 1993; Kollgaard et al. 1996; Rector et 
al. 2000). 

Padovani & Giommi (1995b) have demonstrated that 
the difference in broad-band peaks for HBL and LBL is 
not simply phenomenological. Rather, it represents a fun- 
damental difference between the two sub-classes. In this 
respect, one could expect to find a similar division in the 
FSRQ class - for which, until recently, no evidence ex- 
isted. Indeed, it was suggested by some authors (Sam- 
bruna, Maraschi & Urry 1996), based upon the similari- 
ties of the optical-X-ray broad-band spectral characteris- 
tics of LBL and FSRQ, that no FSRQ with synchrotron 
peak emission in the UV/X-ray band should exist. 

The dichotomy that there existed both high and low 
energy peaked BL Lacs (but virtually no intermediate ob- 
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jects), and only low energy peaked FSRQ was the state of 
the art in our knowledge of blazar spectral energy distri- 
butions (SEDs) in the mid 1990s. Four discoveries have 
drastically changed this picture: 

1. The multi- frequency catalog of Padovani, Giommi, 
& Fiorc (1997), which included all FSRQ known 
before DXRBS and other recent surveys, identified 
more than 50 FSRQ (~ 17% of the FSRQ in their 
catalog) with broad-band spectra similar to those 
of HBL. 

2. About 30% of FSRQ found in the Deep X-ray Radio 
Blazar Survey (DXRBS; Perlman et al. 1998; Landt 
et al. 2001) were found to have X-ray-to-radio lumi- 
nosity ratios, L x /L r , typical of HBL (L x /L r > 10~ 6 
or a rx ^ 0.78), but broad (full width half maximum 
> 2,000 km s _1 ) and luminous (L > 10 43 erg s _1 ) 
emission lines typical of FSRQ. 

3. The discovery of large numbers of "intermediate" 
BL Lacs, objects with broadband properties inter- 
mediate between HBL and LBL, in the DXRBS 
(Perlman et al. 1998; Landt et al. 2001), ROSAT- 
Green Bank (RGB; Laurcnt-Muehleisen et al. 1998) 
and other samples (Nass et al. 1996; Kock et al. 
1996). 

4. BeppoSAX observations have shown that the X- 
ray emission of at least one FSRQ is dominated 
by synchrotron radiation, with a peak frequency 
^peak ~ 2 x 10 16 Hz and steep (a x ~ 1.5) X-ray 
spectrum (Padovani et al. 2002). Two more FSRQ 
show evidence of iy poa k ~ 10 15 Hz. 

The discovery of "X-ray strong" FSRQ (labeled HFSRQ 
by Perlman et al. 1998 to parallel the HBL moniker; the 
"standard" FSRQ would then be the LFSRQ) represents 
a fundamental change in our perception of the broadband 
spectrum of FSRQ, similar to that brought about by the 
X-ray selected BL Lacs discovered in the Einstein Medium 
Sensitivity Survey (EMSS) (Morris et al. 1991; Stocke et 
al. 1991; Rector et al. 2000). A preliminary investigation 
of the properties of HFSRQ in the DXRBS survey was pre- 
sented in Perlman et al. (1998). The discovery of a large 
number of intermediate BL Lacs should have been in many 
ways expected already in the mid 1990s given our knowl- 
edge of the vastly disparate parameter space coverages of 
X-ray and radio surveys (see § 3). 

With these recent discoveries in mind, it is time to 
once again take stock of the broadband spectral proper- 
ties of all blazars. Here we utilize updated identifications 
from DXRBS (most of which were presented by Landt et 
al. 2001), which result in a much larger sample, com- 
bined with a sample of blazars we have extracted from 
RGB (Laurcnt-Muehleisen 1996; Laurcnt-Muchlciscn ct al. 
1998), to attack this problem. A preliminary version of 
some of these findings was presented previously by Perl- 
man (2000). 

The samples arc briefly described in § 2. In § 3 we 
comment on the failure of previous surveys to find X-ray 
strong FSRQ. In § 4 we study the SEDs of our sources and 
discuss the use of the available data to obtain the peak fre- 
quency of synchrotron emission and other characteristics 



of blazars. We discuss the so-called "blazar sequence" sce- 
nario in § 5 and further constraints on blazar emission in 
§ 6, while § 7 summarizes our conclusions. Throughout 
this paper spectral indices are written S v oc v~ a and for 
consistency with previous work the values Hq = 50 km s _1 
Mpc -1 , Om = 0, and f^A = have been adopted. We note 
that our correlations (or lack of) are basically unchanged 
for a cosmological model with Ho = 65 km s _1 Mpc -1 , 
f^M = 0.3, and f^A = 0.7. 

2. THE SAMPLES 

2.1. The DXRBS Blazar Sample 

DXRBS is the result of correlating the ROSAT WGA- 
CAT database (White, Giommi, & Angelini 1995) with 
several publicly available radio catalogs (GB6 and PMN 
at 5 GHz, NORTH20CM at 1.4 GHz), restricting the can- 
didate list to serendipitous flat-spectrum radio sources 
(a r < 0.70). Additionally, a snapshot survey with the Aus- 
tralia Telescope Compact Array (ATCA) was conducted 
for S < 0° to get radio spectral indices unaffected by 
variability (and arc-second positions for the sources south 
of 5 = -40°, the limit of the NRAO-VLA Sky Survey 
[NVSS]; Condon ct al. 1998). The DXRBS X-ray flux 
limits depend on the exposure time and the distance from 
the center of the ROSAT Position Sensitive Proportional 
Counter (PSPC) but vary between - 10~ 14 and - 10~ n 
erg cm -2 s _1 . Radio 5 GHz fluxes reach down to ~ 50 
mJy. 

First results from DXRBS were presented in Perlman 
et al. (1998), while updated identifications were presented 
in Landt et al. (2001). At the time of writing (November 
2002), DXRBS includes 244 blazars: 200 FSRQ (defined 
as broad-line sources with a r < 0.50) and 44 BL Lacs. 
Details on our BL Lac classification are given in Landt et 
al. (2001). Note that Giommi et al. (2002b) have pointed 
out that sources with optical spectrum typical of a galaxy 
but nuclear SED typical of a blazar should probably be 
classified as low-luminosity BL Lacs. Some of our radio 
galaxies might fall in that category. We will address this 
point in future papers. 

We include in this work only sources belonging to the 
complete sample, which fulfills all our selection criteria (see 
details in Landt et al. 2001). This is necessary because 
we want to compare in detail the properties (mainly the 
SED) of DXRBS FSRQ and BL Lacs and that requires 
well-defined samples, not to introduce any bias. Since we 
need the sky coverage to de-convolve the observed distribu- 
tions, we have excluded sources with PSPC center offsets 
in the range 13' — 24', where the sky coverage is difficult to 
determine because of the effects of the spacecraft wobble 
and the rib structure. We are then left with 134 FSRQ 
and 31 BL Lacs, for a total of 165 blazars. Only about 10 
more sources (~ 5% of the total) with a r < 0.50 remain 
to be identified. The current blazar sample is therefore 
highly representative. 

2.2. The RGB Blazar Sample 

Laurent-Muchlciscn ct al. (1997) have described in 
depth the procedures undertaken in constructing the 
ROSAT All Sky Survey (RASS)-Green Bank (RGB) cat- 
alog of radio and X-ray emitting sources. Here we review 
only the essential points. 
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The catalog was constructed by correlating the RASS 
with a radio catalog based on the 1987 Green Bank 5 
GHz survey maps (Gregory & Condon 1991; Gregory et 
al. 1996). VLA observations of all the 2,127 radio/X-ray 
matches within 100" were obtained to derive arc-second 
positions (Laurent-Muehleisen et al. 1997). For the 1,567 
sources with radio/X-ray offset < 40" optical counter- 
parts were identified using the Automatic Plate Measur- 
ing (APM) scans (Irwin, Maddox, & McMahon 1994), 
which give O (blue) and E (red) magnitudes. Laurent- 
Muehleisen et al. (1997, 1998) give many new object 
identifications based upon spectra taken at various opti- 
cal telescopes. The identifications published by Laurent- 
Muehleisen et al. (1997, 1998) include a comprehensive list 
of BL Lacs, as well as many radio-loud, broad-line objects 
(i.e., radio- loud quasars); however, since radio spectral in- 
dex was not included in the search criteria for RGB no list 
of FSRQ was compiled there. 

Some indication that X-ray strong FSRQ were being 
found in the RGB sample could be seen from the distri- 
bution of X-ray to radio flux ratios of the newly identified 
objects (Laurent-Muehleisen et al. 1998). We therefore de- 
cided to use the information in the RGB along with other 
public surveys to create a list of all blazars (BL Lacs and 
FSRQ) in the full RGB sample. We used Tables 3.3 and 
3.4 of Laurent-Muehleisen (1996) and Table 3 of Laurent- 
Muehleisen et al. (1998) to search for both BL Lacs and 
FSRQ as follows. First, the RGB sample was cross- 
correlated with the AGN catalog of Padovani, Giommi, 
& Fiore (1997) to find previously known blazars. Then, 
we used the identifications given by Laurent-Muehleisen 
et al. (1998) to include newly identified sources (note that 
these are limited to objects with O < 18.5). This resulted 
in a total of more than 600 classified RGB sources. Being 
based on the RASS, the RGB covers a very different X-ray 
flux range: X-ray fluxes range between ~ 2 x 10~ 13 and 
~ 10~ 10 erg cm~ 2 s _1 . Radio 5 GHz fluxes reach down to 
~ 25 mJy, roughly a factor of 2 fainter than DXRBS at 
similar declinations. 

Radio spectral indices, necessary to determine if a 
broad-line radio source is an FSRQ or not, were derived by 
cross-correlating the RGB with the NVSS. For consistency 
with DXRBS, we used 5 GHz fluxes from the GB6 survey 
(and not from the dedicated VLA observations) for the 
RGB sources. In doing the cross-correlation, then, the 1.4 
GHz flux from all NVSS sources within a 3' radius (cor- 
responding roughly to the beam size of the GB6 survey) 
was summed up (we have already described the rationale 
connected with doing this in Pcrlman et al. 1998). We 
then isolated a sample of FSRQ (a r < 0.5). 

We checked that the spectral indices derived from the 
GB6/NVSS radio data were reliable in the following way. 
Three hundred and forty-five radio-loud RGB sources had 
a value of the radio spectral index in the AGN catalog 
of Padovani, Giommi, & Fiore (1997), generally obtained 
from non-simultaneous, single-dish measurements around 
a few GHz. There is a very strong (P > 99.99%) linear 
correlation with a slope ~ 0.9 between the two spectral 
indices, with a mean difference Aa = 0.08 over a range 
of ~ 2.5. This shows that the method we used to esti- 
mate the radio spectral index for the RGB sample is quite 
robust. 



The RGB blazar sample thus includes 362 blazars: 233 
FSRQ and 129 BL Lacs. To our knowledge the RGB sam- 
ple is <~ 44% identified. This fraction goes up to ~ 76% for 
the part of the sample with O < 18.5. Considering only 
the sources with a r < 0.5, these fractions do not change 
much, being 49% and 73% respectively. 



2.3. The Joint Sample 

Overall, the two samples used in this paper include 497 
distinct blazars, 342 FSRQ and 155 BL Lacs (30 objects 
are in common). In terms of its range of properties, size, 
depth, and selection criteria this ensemble of objects rep- 
resents a unique sample with which to address some of the 
open questions of blazar research. 

As an initial step towards studying the broad-band 
properties of our sources, we first derive their a ox , a ro , and 
a rx values. These are the usual rest-frame effective spec- 
tral indices defined between 5 GHz, 5,000 A, and 1 keV. 
X-ray and optical fluxes have been corrected for Galac- 
tic absorption. The effective spectral indices have been 
k-corrected using the appropriate radio, optical, and X- 
ray spectral indices. X-ray spectral indices are available 
for all DXRBS sources from hardness ratios, as described 
in Landt et al. (2001). For the RGB objects we obtained 
X-ray spectral indices from hardness ratios for 73% and 
89% of FSRQ and BL Lacs respectively from Brinkmann 
et al. (1997), Laurent-Muehleisen et al. (1999), Reich et al. 
(2000), and WGACAT. In the remaining cases we assumed 
a x = 1.2 following Laurent-Muehleisen et al. (1998). Val- 
ues for the optical spectral indices were assumed to be the 
typical ones for the various sub-classes. Finally, optical 
magnitudes for broad-line sources were also corrected for 
the presence of emission lines according to the prescription 
of Natali et al. (1998). 

It is worth mentioning that, unlike DXRBS, all new 
RGB identifications have O < 18.5. This has important 
selection biases associated with it (§§ 3,4). In particular, 
an optical magnitude limit has the effect of imposing an 
artificial boundary on the regions of parameter space to 
which a survey is sensitive. Namely, only sources with 
a m < a ro (lim) and a ox > a ox (lim), where these limiting 
values depend on radio and X-ray flux, will be included. 
At the survey radio limit, for example, only sources with 
a ro ^5 0.4 will be included. We will discuss the effect of 
this limitation in § 4. 



3. THE DISCOVERY OF X-RAY STRONG FSRQ 

The findings of Padovani, Giommi, & Fiore (1997), Perl- 
man et al. (1998), and Padovani et al. (2002), make it clear 
that detecting X-ray strong FSRQ requires large and/or 
deep samples. It was already known in 1995 that radio- 
selected samples find more frequently objects with broad- 
band spectra peaking in the IR/optical. The oft-repeated 
separation of HBL and LBL on the a ox , a ro plane (e.g., 
Padovani & Giommi 1995b), is a direct by-product of these 
disjoint survey methods. 
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Fig. 1.— The a TO , a ox plane for the DXRBS (top) and RGB (bot- 
tom) samples. Effective spectral indices are defined in the usual way 
and calculated between the rest-frame frequencies of 5 GHz, 5000 
A, and 1 kcV. Filled circles represent the DXRBS and RGB FSRQ, 
open squares represent the 1 Jy FSRQ, while crosses are the DXRBS 
and RGB BL Lacs. The dashed lines represent, from top to bottom, 
the loci of o rx = 0.85, typical of 1 Jy FSRQ and LBL, a rx = 0.78, 
the dividing line between HBL and LBL, and c? rx = 0.70, typical of 
RGB BL Lacs. The regions in the plane within 2cr from the mean 
ctro , c*o X , and o rx values of LBL and HBL are indicated by the solid 
lines and marked accordingly. 

With this observation in mind, one might expect a sim- 
ilar selection effect to have been present for FSRQ. Yet 
up until the mid-1990s, no X-ray survey had looked for 
FSRQ, even though re-examinations of the Einstein EMSS 
and Slew Survey databases have revealed samples of FSRQ 
in both (Perlman 2000; Perlman et al. 2001; Wolter & 
Celotti 2001). The only flux-limited samples of FSRQ 
which existed prior to that time, in fact, were "classical" 
high flux limit, radio-selected samples, such as the 1 Jy, 
produced by surveys which also contained few, if any HBL 
(e.g., the 1 Jy sample has only 2 HBL out of 34 BL Lacs 
total, Padovani & Giommi 1995b). Thus it was perhaps 
not a surprise that these radio-selected samples contained 
very few X-ray strong FSRQ. 

Figure 1 shows the distribution of our sources in the 
a x,a ro plane. The figure shows also the 1 Jy FSRQ, 
which occupy a region of a ox ,a ro parameter space with 
a rx similar to that typical of LBL (marked in the figure). 
FSRQ with low a rx (< 0.78, roughly equivalent to the 
HBL/LBL division (although see § 4)', or L x /L r > 10~ 6 ) 
constitute only ~ 5% of the 1 Jy sources with X-ray data. 
Importantly, none of the 1 Jy FSRQ fall in the region of 
the plane within 2a from the mean a ro , a ox , and a rx val- 
ues of HBL, the "HBL box" , derived by using all HBL in 
the multi-frequency AGN catalog of Padovani, Giommi, & 
Fiore (1997) 5 . The fainter radio/X-ray selected DXRBS 
and RGB FSRQ, on the other hand, reach much lower 



values of a rx and many sources "invade" the HBL region. 
Indeed, as can be seen by comparing the two panels of 
Figure 1, there is a progression of a ox ,a ro from 1 Jy to 
DXRBS to RGB, with the exception of the half of the 
HBL box to the left of the diagonal line that extends from 
(0.8,0.55) to (1.4,0.2). We will return to this subject in 
§4. 
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Fig. 2. — The effective radio-optical spectral index a IO vs. the 5 
GHz radio flux for about 900 FSRQ from the multi-frequency cata- 
log of Padovani, Giommi, & Fiore (1997) (crosses) and the DXRBS 
FSRQ discussed in this paper (filled points). 

Another way of looking at selection effects is to consider 
Fig. 2, which plots a m vs. the 5 GHz radio flux for about 
900 FSRQ from the comprehensive multi-frequency cata- 
log of Padovani, Giommi, & Fiore (1997). The observed 
trend between radio flux and optical-radio spectral index, 
with stronger radio sources having steeper a ro values, is 
likely due to two separate effects: 1. the relatively bright 
optical flux limits of the various samples which were in- 
cluded in the catalog, responsible for the lack of sources in 
the upper left part of the diagram. Indeed DXRBS FSRQ, 
which reach V > 23, start to fill this region (filled points 
in Fig. 2); 2. the fact that well-known, "classical" FSRQ, 
all relatively strong radio sources (3C 273 and 3C 279, for 
example, have 5 GHz radio fluxes ~ 16 and 40 Jy respec- 
tively) do have steep a ro , which accounts for the lack of 
sources in the lower right part of the diagram. Since HF- 
SRQ need to have a ro < 0.5 (see Fig. 1) they will start to 
be revealed only at fainter (< 200 — 300 mJy) radio fluxes, 
and it is only by conducting relatively faint radio surveys 
with X-ray information, like ours, that these sources can 
be identified in reasonable numbers. 

There could be two possible reasons for the observations 
pointed out in the previous two paragraphs. Either some 
selection effect is at work, perhaps induced by a correlation 



5 X-ray data are available for ~ 65% of 1 Jy FSRQ. The sources without X-ray data have c* ro > 0.63 and therefore are all outside of the "HBL 
box". 
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between the fraction of HFSRQ and flux, with HFSRQ be- 
ing intrinsically more numerous at fainter fluxes/powers, 
or these sources are intrinsically rare and their detection 
in sizeable numbers requires large samples. We will ad- 
dress this issue in future DXRBS papers. Note that simi- 
lar questions are perfectly appropriate as well for BL Lacs 
and these are the subject of some debate (see, for example, 
Fossati 2001, Giommi et al. 2002c). 



4. SEDS AND THE NATURE OF BLAZAR EMISSION 

In the previous section we have shown that recent, rel- 
atively deep, blazar samples with both radio and X-ray 
information have detected the hitherto unknown class of 
FSRQ with effective spectral indices typical of HBL. Based 
on our knowledge of BL Lacs, we expect these sources to 
have a synchrotron peak frequency t'peak in the UV/X-ray 
band. Similarly, we would expect that the "intermedi- 
ate" BL Lacs in our sample should have synchrotron peak 
frequencies intermediate between those seen for HBL and 
LBL, as demonstrated for S5 0716+714 by Giommi et al. 
(1999). We then need to study the SEDs of our sources. 



4.1. a rx Distributions 

The HBL/LBL division is generally done in terms of the 
X-ray-to-radio flux ratio or effective radio-X-ray spectral 
index a IX (Padovani & Giommi 1995b). While this pa- 
rameter is broadly related to z/ pca k (Padovani & Giommi 
1996; Fossati et al. 1998; but see § 6.2), it is much easier 
to derive. Fig. 1 shows that the DXRBS FSRQ and BL 
Lacs seem to cluster at the edge of the HBL region, while 
the RGB FSRQ and BL Lacs are moving in progressively. 
This is due at least partly to the different flux limits of the 
two surveys, which sample different regions of parameter 
space. It does appear, however, that while the RGB BL 
Lacs populate a much larger area of the HBL region than 
DXRBS BL Lacs, this is not the case for FSRQ. Inclusion 
of the unclassified RGB sources with optical information 
does not alter this conclusion (see also below). Optically 
fainter sources will have relatively flat a ox but most impor- 
tantly relatively steep a ro and so will be out of the HBL 
box. 

Fig. 3 shows the fractional ct rx distribution of all 
DXRBS and RGB blazars. It is important to note that due 
to the serendipitous nature of DXRBS (and other, previ- 
ous serendipitous surveys, including the EMSS), the area 
in which faint X-ray sources could be detected is smaller 
than that for brighter X-ray sources. Thus, in order to 
compare directly distributions of parameters related to X- 
ray flux, we need to "correct" by using the appropriate sky 
coverage - a necessary step which is seldom done - oth- 
erwise faint sources would be under-represented. We have 
already done this for the redshift distribution in Landt et 
al. (2001) and we refer the reader to that paper for fur- 
ther details. a rx is obviously related to X-ray flux so we 
have de-convolved its distribution, shown in Fig. 3. We 
follow this practice for all comparisons involving DXRBS 
and any other serendipitous survey. 
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Fig. 3. — The fractional distribution of the X-ray-radio spectral 
index ct rx for FSRQ (solid lines) and BL Lacs (dashed lines) for the 
DXRBS (top) and RGB (bottom) samples. The dotted line repre- 
sents the Q rx distribution of RGB FSRQ and unclassified radio-loud 
sources with o r < 0.5. The DXRBS distributions have been de- 
convolved with the appropriate sky coverage. Error bars represent 
the la range based on Poisson statistics. See text for details. 



A few points can be inferred from Fig. 3. The (a rx ) 
values for DXRBS FSRQ and BL Lacs are similar, ~ 0.83, 
but the variances are different at the 99.96% level (Stu- 
dent's t-test). Moreover, the two distributions are differ- 
ent at the 99.4% level according to a Kolmogorov-Smirnov 
(KS) test. The BL Lac a TX distribution, in fact, is broader 
than that of FSRQ, with a peak at a rx ~ 0.9 and a tail 
reaching relatively low values (~ 0.55). The skewness of 
the BL Lac distribution is negative and equal to ~ 3.4 
times its standard deviation, unlike the FSRQ distribu- 
tion for which the skewness is consistent with being zero 
(Press et al. 1986). The different nature of the BL Lac 
and FSRQ a rx distributions is confirmed by the use of 
the kmm algorithm, developed by Ashman, Bird, & Zcpf 
(1994), which computes for a given distribution the confi- 
dence level at which a single Gaussian fit can be rejected 
in favor of a double Gaussian fit. While for the DXRBS 
FSRQ we find that a single Gaussian provides a good fit 
to the a IX distribution, this possibility is rejected at the 
99.9% level for the BL Lacs. The - 10 still unidentified 
sources with a r < 0.5 are X-ray faint and therefore have 
relatively large (a rx ) ~ 0.9. Based on our identification 
statistics, we expect most of them to be FSRQ, thereby 
somewhat reducing the difference in the a IX distributions 
at large values in Fig. 3. The fraction of sources with 
L x /L r > 10~ 6 , i.e., a IX < 0.78 is - 28% and - 25% for 
the BL Lacs and FSRQ respectively. More interesting is 
the fraction of sources which fall into the region typical of 
HBL, the "HBL box". This is - 15% and - 9% for the 
BL Lacs and FSRQ respectively. For RGB the a IX distri- 
butions for FSRQ and BL Lacs are also significantly dif- 
ferent (P > 99.99%), but in this case the mean values are 
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also different, with (a rx ) ~ 0.80 and ~ 0.68 respectively. 
The fraction of sources with L x /L r > 10~ 6 is ~ 76% and 
~ 35% for the BL Lacs and FSRQ respectively. The frac- 
tion of sources falling into the "HBL box" is ~ 60% and 
~ 27% for the BL Lacs and FSRQ respectively. While 
RGB BL Lacs reach more extreme (i.e., much flatter) a rx 
values than DXRBS ones, that is not the case for FSRQ. 

Could this difference between RGB FSRQ and BL Lacs 
be explained as a selection effect? A large fraction (<~ 60%) 
of the most X-ray loud RGB BL Lacs (a IX < 0.7) were al- 
ready known because of dedicated surveys (EMSS, Slew, 
HEAO-1). Given the fact that the identification fraction of 
the flat-spectrum sources in RGB is quite low (~ 49%) and 
that the identified sources in RGB have mostly O < 18.5, 
it could be that the most extreme RGB FSRQ are still 
unclassified. We have then included with the RGB FSRQ 
the ~ 360 unclassified radio-loud sources with a r < 0.5 
(assuming z = 1 for the K-correction) . Note that this is 
an extreme assumption as a substantial fraction of these 
could be BL Lacs (their current fraction in the a r < 0.5 
sample is ~ 30%). Although the a IX distribution now 
shifts to lower values, with (a IX ) ~ 0.74, it is still signif- 
icantly different (P > 99.99%) from that of the RGB BL 
Lacs. The fraction of sources with L x /L r > 10~ 6 rises 
from ~ 35% to ~ 62%, while that which falls into the 
"HBL box" (amongst the sources with optical informa- 
tion; ~ 40% of the unclassified sources are empty fields) 
rises from ~ 27% to ~ 38%. As discussed above, this 
fraction should probably be considered an upper limit. 

We point out that the DXRBS and RGB a IX distribu- 
tions for both FSRQ and BL Lacs peak at ~ 0.7 — 0.8, 
that is around the HBL/LBL dividing line. The previously 
suggested dichotomy in the BL Lac class then, with most 
sources at the two ends of the distribution, was simply 
due to the selection effects induced by combining widely 
disconnected samples selected in different bands. 

We can check what happens at even larger X-ray-to- 
radio flux ratios by using two other samples: the EMSS 
and the Sedentary survey. The EMSS (e.g., Maccacaro 
et al. 1994) is an X-ray selected sample of sources dis- 
covered in 1,435 Einstein Imaging Proportional Counter 
(IPC) fields. It reaches ~ 5 x 10 -14 erg cm -2 s _1 in the 
0.3 — 3.5 keV band but the area covered is a strong func- 
tion of X-ray flux (Gioia et al. 1990). The sample is largely 
identified (completely so down to ~ 2 x 10 -13 erg cm~ 2 
s _1 ) and comprehensive, dedicated radio observations pro- 
vide detections down to ~ 1 mJy levels so that quite high 
values of X-ray-to-radio flux ratios can be reached. 

Two papers have constructed samples of FSRQ within 
the EMSS. Wolter & Celotti (2001) constructed a sample 
of EMSS radio-loud quasars by using a cut a IO > 0.35. 
We have used the values for the thirteen FSRQ (a r < 0.5) 
given in their paper. Perlman et al. (2001) did a simi- 
lar search which yielded those objects plus two more with 
0.2 < a ro < 0.35 (for consistency with the commonly used 
definition of radio-loud sources; these two sources happen 
to have the two lowest a IX values). Using these objects, 
we have constructed the a IX distribution of these EMSS 
FSRQ, again taking into account the effect of the sky cov- 
erage. We have also done the same for the EMSS BL 
Lac sample, using the revised list and data of Rector et 
al. (2000). The results are as follows: the a rx distribu- 



tions for FSRQ and BL Lacs are still significantly different 
(P > 99.99%), with mean values (a rx ) ~ 0.71 and ~ 0.60 
respectively. The minimum value reached by EMSS FSRQ 
is a rx ~ 0.63, while that of EMSS BL Lacs is ~ 0.49. The 
fraction of sources with L x /L r > 10~ 6 , once the effect of 
the sky coverage is taken into account, is 100% and ~ 78% 
for the BL Lacs and FSRQ respectively. The fraction of 
sources which fall into the "HBL box" is the same, namely 
100% and - 79% for the BL Lacs and FSRQ respectively. 
So, although the a rx gap between the two classes gets nar- 
rower, the fact remains that FSRQ do not reach the same 
values as BL Lacs. 

An even more extreme sample is the Sedentary survey 
(Giommi, Menna, & Padovani 1999), an X-ray/radio se- 
lected sample based on the RASS Bright Source Catalog 
(RASSBSC) and the NVSS, which has a cut at a rx < 0.56. 
The current sample is ~ 90% identified so our conclu- 
sions should be relatively stable. The number of radio- 
loud broad-lined sources in the sample is currently around 
20 (or ~ 12%), but most of these sources are very close 
to the radio-loud/radio-quiet dividing line, both in terms 
of their a ro values and their radio powers. Some of these 
sources could therefore have their X-ray emission domi- 
nated by thermal processes and would not be the coun- 
terparts of HBL. We will discuss the presence of HFSRQ 
in the Sedentary survey in a future paper (Giommi et al., 
in preparation). Table 1 summarizes the results of these 
comparisons, in terms of mean a rx values and fractions of 
HBL and HFSRQ in the DXRBS, RGB, and EMSS sam- 
ples. 

Our results reaffirm the fact that the percentage of 
HBL/HFSRQ in a sample is heavily dependent on the sur- 
vey flux limits or, alternatively, its position on the X-ray- 
radio flux plane (see Figs. 1 and 2 of Padovani (2002)). 
Tab. 1 shows very explicitly how these percentages change 
by moving from DXRBS, a survey whose limits are in the 
"L" zone, to RGB and EMSS, whose limits move progres- 
sively deeper into the "H" zone. 

4.2. Synchrotron Peak Frequency 

To study in more detail the synchrotron peak frequen- 
cies of our sources we have used the multi-frequency infor- 
mation at our disposal to extract non-simultaneous SEDs 
for our blazars. We restrict ourselves to the DXRBS sam- 
ple as the determination of ^ poa k requires a lot more effort 
than a rx and we did not think it was worth deriving for 
RGB given its large incompleteness. 

We determined f poa k for all 165 DXRBS blazars by 
applying an homogeneous synchrotron - inverse self- 
Compton (SSC) model in the log v - log vf v plane to 
the following data: 

1. radio fluxes at two different frequencies (non- 
simultaneous 1.4 and 5 GHz data for objects with 
S > —40°; simultaneous 5 and 8.6 GHz data for 
objects with S < -40°; see Landt et al. 2001); 

2. optical flux, in the V band, derived from Palo- 
mar Observatory Sky Survey (POSS I) O and/or 
E magnitudes as described in Giommi, Menna, & 
Padovani (1999); fluxes were corrected for Galactic 
absorption and for the presence of emission lines 
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Table 1 
Q rx Statistics. 
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RGB 


0.803 ± 0.005(> 0.74) 
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27%(< 38%) 


60% 
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0.71 ±0.02 


15 


0.60 ± 0.01 


41 


78% 


100% 


79% ~ 
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(for FSRQ) according to the prescription of Natali 
et al. (1998); 

3. unabsorbed ROSAT X-ray flux at 1 keV, derived 
from the broad-band flux and the X-ray spectral 
index a x ; 

4. Two Micron All Sky Survey (2MASS) data (if avail- 
able; Cutri ct al. 2000); 

5. any other NASA/IPAC Extragalactic Database 
(NED) data. 

We note that here and in the following we deal with the 
rest-frame peak frequency, equal to i/ pea k(obs) x (1 + z). A 
rcdshift of 0.4 was assumed for the BL Lacs without one. 

While we have optical spectra for all objects identified in 
Pcrlman et al. (1998) and Landt et al. (2001), we did not 
use those spectra to estimate peak frequencies for two rea- 
sons. First, the DXRBS includes also a significant number 
of previously known sources (about 40% for the sample 
used in this paper), for which optical spectra are not nec- 
essarily available in digital format. We felt it would be 
inconsistent to use optical spectra for some objects but 
not for others. Second, for many objects the spectra were 
taken with a slit that was not at parallactic angle due to 
the instrumental setups used (see Perlman et al. 1998 and 
Landt et al. 2001). This results in a flux loss due to atmo- 
spheric differential refraction which is more pronounced in 
the blue part of the spectrum and would therefore signifi- 
cantly affect the optical/UV slope. 

The SSC model was been adapted from Tavecchio, 
Maraschi & Ghiscllini (1998) and assumes that radiation 
is produced by a population of relativistic electrons emit- 
ting synchrotron radiation in a single zone of a jet that is 
moving at relativistic speed and at a small angle to the 
line of sight. These photons are subsequently scattered 
by the same electrons to higher energies via the inverse 
Compton process (note that we ignore Comptonization of 
external photons, which only affects the SED at 7-ray ener- 
gies; e.g., Ghisellini et al. 1998). The physical parameters 
that define the model are the jet radius, the Doppler fac- 
tor 6, the magnetic field B, and four spectral parameters 
of the electron population, assumed to follow a power-law 
distribution which breaks sharply above a given energy: 
the normalization, the two spectral slopes, and the break 
energy. The Klein-Nishina cross section is used in the com- 
putation of the Compton scattering. 

Our aim is to derive a best estimate of f poa k- Given 
that the SED of most of our sources is not well sampled 
we cannot fully constrain the model parameters. However, 
the available range was limited by physical considerations, 
namely: 6 ~ 10 - 20, B ~ 0.1 - 5 G, jet radius « 10~ 3 



parsec, and a low-energy slope of the electron distribution 
which produces a relatively flat (< 0.7) spectral slope in 
the radio band, consistent with the sample definition. We 
have also been guided by the work done by some of us 
using the same model to derive f pca k for a large sample 
of bright blazars observed by BeppoSAX (Giommi et al. 
2002a) and from the NVSS-RASS 1 Jy survey (Giommi et 
al. 2002b). Both samples, in fact, include many "classi- 
cal" blazars with well sampled SEDs for which the model 
parameters were well constrained. In any case, our experi- 
ence shows that ^ pC akj a combination of B, 8, and electron 
break energy, is relatively well constrained even for rel- 
atively large variations of these parameters, which affect 
much more the high-energy/inverse Compton part of the 
SED (e.g., Massaro et al. 2002). 

FSRQ can have a disk (thermal) component in the 
optical/UV band but this, on average, makes up only 
~ 15% of the continuum emission in DXRBS FSRQ 
(D'Elia, Padovani, & Landt 2003), and therefore should 
not strongly affect our derivation of !/ pC akj at least in a 
statistical sense, although there are individual exceptions. 
This is also consistent with the relatively steep optical/UV 
slopes of these objects (a ~ 1.2; Landt et al., in prepara- 
tion). We think it is also unlikely that the thermal com- 
ponent can make a substantial contribution in the X-ray 
band, as our FSRQ have typical redshifts ~ 1.5, which 
implies that the ROSAT band corresponds to 0.25 — 6 keV 
rest-frame. By comparison, there is very little evidence 
for thermal optical/UV emission due to an accretion disk 
in BL Lacs, although Corbett et al. (2000) have suggested 
that accretion disk illumination is responsible for the vari- 
ations in Ha luminosity in BL Lacertae itself. In any case, 
however, the line luminosity of DXRBS BL Lacs is much 
lower than that of FSRQ (Landt et al., in preparation), 
so that any contribution to the total optical/UV flux is 
negligible for our purposes. 

Previous works have derived f pea k by fitting analytical 
functions, such as a parabola or a third-degree polynomial, 
to the SED of blazars (e.g., Sambruna et al. 1996; Fossati 
et al. 1998). Although the latter alternative would also 
allow us to take into account the spectral upturn which 
is present when the X-ray flux is dominated by inverse 
Compton emission, we chose instead to derive f poa k using 
a physical model (SSC) which is widely accepted as being 
responsible for the broad-band emission of blazars. We 
believe that our approach, although more complex and 
time consuming, is more robust especially when dealing 
with sparsely sampled SEDs, as we are guided by physics 
rather than just analytical fitting. We also note that our 
multifrcquency data are not simultaneous and therefore 
our fits are also affected by variability. 

As examples, in Figs. 4 and 5 we show some representa- 
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Fig. 4. — Representative synchrotron self-Compton fits to the spectral energy distributions of DXRBS FSRQ. Filled points denote DXRBS 
data, empty points denote data from NED and 2MASS. The synchrotron peak frequency for the two sources are t'peak = 5.5 X 10 15 Hz (WGA 
J0447.9-0322, z = 0.774, left) and ^ pcak = 1.7 x 10 13 Hz (WGA J0510.0+1800, z = 0.416, right). Sec text for details. 




Fig. 5. — Representative synchrotron self-Compton fits to the spectral energy distributions of DXRBS BL Lacs. Filled points denote 
DXRBS data. The synchrotron peak frequency for the two sources are u pcak = 1.1 X 10 15 Hz (WGA J0043.3-2638, z = 1.002, left) and 
^pcak = 2.3 x 10 13 Hz (WGA J0533.6-4632, z = 0.332, right). See text for details. 



tive fits to the SEDs of DXRBS FSRQ and BL Lacs, span- 
ning a wide range of ^ pca k- We stress again that the sam- 
pling of the SEDs is scanty, being based in many cases only 
on a few, non-contemporaneous data points, so that the 
precise value of f poa k for a given source can be somewhat 
uncertain. However, we do believe that its statistical use is 
warranted. This is confirmed by plotting the optical/UV 
slopes, which are measured for 68 (~ 51%) of our FSRQ 
and did not enter in our derivation of the synchrotron 
peak frequency, vs. ^ pca k- A very strong anti-correlation 
(P = 99.96%) is observed between the two quantities, as 
expected if our ^ pca k estimates were indeed reliable. For 
relatively small f pca k values, in fact, one expects the opti- 



cal band to sample the steep tail of the synchrotron com- 
ponent, while when the synchrotron peak moves to pro- 
gressively larger energies the optical slope becomes flatter 
as the optical band samples the lower frequency, flatter 
parts of the synchrotron emission. 

Figure 6, which shows the ^ pca k distribution of all our 
sources, illustrates the fact that indeed, the SEDs of some 
DXRBS FSRQ peak in the UV/X-ray band. As in the case 
of a rx , the f poa k distributions are also corrected for the ef- 
fect of the sky coverage. The FSRQ distribution ranges 
between 10 12 and 10 16 Hz, with (z/ pcak ) = io 13 - 70±0 06 Hz, 
while the BL Lac distribution ranges between 6 x 10 12 
and 8 x 10 16 Hz, with (j/ pcak ) = 10 141±01 Hz. On aver- 
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age, then, BL Lacs have a synchrotron peak frequency a 
factor ~ 2.5 larger than that of FSRQ. The two distribu- 
tions are also different at the 99.97% level according to a 
KS test. As was the case for a rx , the BL Lac distribution 
is broader than the FSRQ one, with a tail reaching higher 
values. For example, while ~ 64% and ~ 8% of BL Lacs 
have J^pcak > 10 14 and 10 15 Hz, these fractions go down to 
~ 31% and - 4% for FSRQ. 

We cannot think of any reason why our method should 
result in lower ^ pca k for FSRQ. Quite the opposite: both 
the thermal component (which is small, as discussed 
above) and the emission lines contamination of the op- 
tical flux (for which we correct anyway) would move the 
peak to higher values (at least with relatively low t'peak)- 
We therefore regard this difference as real. 

It is likely that this difference in z/ pca k, actually, is even 
larger than shown here. Due to its serendipitous nature, 
DXRBS is subject to incompleteness at high fluxes, since 
bright targets are excluded by definition not to bias the 
sample. We can quantify this in the following way. Ap- 
proximately 50 BL Lacs with radio flux above our thresh- 
old, mostly from the 1 Jy, Slew, and EMSS samples, and 
off the Galactic plane have been observed by ROSAT as 
targets. This corresponds to a surface density «2x 1CP 3 
deg -2 , i.e. roughly 3 — 4 sources which could not be in- 
cluded in our survey. This effect is more severe at high 
^peak values where only a few objects have been detected 
(see Fig. 6) and the loss of even a small number of sources 
could make a major difference. Note in fact that many 
of these BL Lacs are HBL. As regards FSRQ, although 
the numbers are somewhat larger, the effect at high f pea k 
values is zero because none of the previously known (and 
therefore likely ROSAT targets) FSRQ had high ^ poa k. 
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Fig. 6. — The distribution of the synchrotron peak frequency for 
FSRQ (solid line) and BL Lacs (dashed line) for the DXRBS sam- 
ple. The distributions have been de-convolved with the appropriate 
sky coverage. Error bars represent the la range based on Poisson 
statistics. See text for details. 
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Fig. 7. — Radio power at 5 GHz vs. the synchrotron peak 
frequency for FSRQ (filled points) and BL Lacs (crosses) for the 
DXRBS sample. The dotted lines denote the two quadrants (top- 
left and bottom-right) occupied by the sources studied by Fossati et 
al. (1998). See text for details. 



While already in DXRBS then there is evidence that 
FSRQ do not appear to reach the more extreme a rx and 
^pe a k values of BL Lacs, this tendency is confirmed by 
studying the RGB and EMSS samples, which sample the 
blazar population deeper into the region of lower a rx . The 
fact that the Sedentary survey, which samples the extreme 
HBL zone, has found very few, if any, FSRQ, is also sug- 
gesting that there might be an intrinsic limit to the ^ poa k 
values which can be reached by strong-lined blazars. 

5. THE BLAZAR SEQUENCE 

Fossati et al. (1998) and Ghiscllini et al. (1998) have pro- 
posed that some blazar properties can be accounted for by 
an inverse correlation between intrinsic power and the syn- 
chrotron peak frequency, the so-called "blazar sequence" . 
The peak of the emission is related to the electron energy, 
as ^ P c a k oc -B7b rea k, c witn 7brc a k,c a characteristic electron 
energy which is determined by a competition between ac- 
celeration and cooling processes. Therefore, less powerful 
sources (where the energy densities are relatively small) 
should reach a balance between cooling and acceleration 
at larger f poa k, while in more powerful sources there is 
more cooling and the balance is reached at smaller ^ pca k- 

Figure 7 plots radio power at 5 GHz versus ^ pea k for the 
DXRBS sources. The dotted lines denote the two quad- 
rants (top-left and bottom-right) occupied by the sources 
studied by Fossati et al. (1998), which belonged to the 1 
Jy and Slew BL Lac samples and the 2 Jy FSRQ sample. 

A few points can be made about this figure. First, as 
already shown in Fig. 6, DXRBS BL Lacs reach f pea k 
values slightly higher than DXRBS FSRQ. (Based on the 
analysis shown in § 4.1 this difference is likely to become 
more pronounced for more X-ray extreme surveys such as 
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EMSS.) Second, DXRBS sources are starting to occupy re- 
gions of this plot (top-right and particularly bottom-left) 
which were "empty" in the original plot of Fossati et al. 
(1998). DXRBS reaches lower radio powers than the 1 Jy 
sample and therefore detects low-power LBL. In particu- 
lar, out of the 21 BL Lacs with ^ pca k < 10 15 5 Hz and 
redshift information, 7 (or 33%) "invade" the low-power 
part (L r < 10 25 - 3 W/Hz) of the plot. Third, no correla- 
tion (P ~ 93%) is present between radio power and v pe &k 
for the whole sample, or for the FSRQ and BL Lac sam- 
ples separately. We also note that the scatter in the plot is 
very large, reaching four orders of magnitude in power for 
10 13 < ^p ak S 10 15 Hz. Therefore, by using an homoge- 
neous, well-defined sample which includes both FSRQ and 
BL Lacs over a relatively wide region of parameter space, 
the Fossati et al. correlation is not confirmed. Fourth, 
an upper envelope, however, seems to be present in the 
right part of the diagram. For example, all sources with 
L r > 10 27 - 5 W Hz -1 have f poak < 10 15 Hz while sources 
with L r > 10 28 W Hz" 1 have i/ poak < 10 14 Hz. 

We have checked the radio power and ^ pca k values and 
distributions for DXRBS blazars in and out of the HBL 
box. We find that the two classes have indistinguishable 
radio powers but significantly different synchrotron peak 
frequencies, with mean values (^ pca k) ~ 10 15 ' 2±f)1 and 
~ io 13 - 66±0 05 Hz respectively, again in contrast to the 
proposed correlation. 

We note that the lack of high z/ pca k high radio power 
blazars could also be due to a selection effect, as discussed 
by Giommi et al. (2002c) (see also (Padovani et al. 2002)). 
This combination, in fact, implies a dominance of non- 
thermal emission over the host galaxy and emission line 
components, making the redshift determination very hard 
if not outright impossible. Objects of this kind, which 
would populate the upper right part of Fig. 7, would then 
be excluded. We note, however, that only four DXRBS 
BL Lacs have no redshift and ^ poa k > 10 14 5 Hz. 




It is interesting to explore possible correlations between 
^pe a k and other powers. The correlation suggested by Fos- 
sati et al. (1998) and Ghiscllini et al. (1998) was between 
synchrotron peak frequency and intrinsic power. As our 
sources are all flat-spectrum their radio power is strongly 
affected by beaming and this could influence the inter- 
pretation of Fig. 7. We have then evaluated two intrin- 
sic powers for our FSRQ, namely the Broad Line Region 
(BLR) luminosity, Lblr, following Celotti, Padovani, & 
Ghiscllini (1997), and the kinetic jet power, Lj Ct , following 
D'Elia, Padovani, & Landt (2003). 

Lblr is an isotropic quantity, related to the ionizing, 
disk (thermal) emission via the covering factor / cnv , i.e., 
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Padovani, & Landt 2003). Figure 8 plots Lblr vs. v pe a.k 
for the 94 (~ 70%) DXRBS FSRQ for which wc could find 
the relevant BLR line fluxes. We do not include here BL 
Lacs because only 6 of them display the lines needed to 
derive Lblr and even in those cases the lines are narrow. 
No hint of an inverse correlation between Lblr, propor- 
tional to disk power, and ^ poa k is present in the figure, 
the opposite is rather apparent. There is in fact a strong 
correlation (P > 99.99%), with a large scatter, between 
the two quantities, with Lblr °c ^ P o a k ±0 U - ^ could be 
argued that the stronger the disk emission, the larger its 
contribution to the optical/UV flux, and the higher the 
estimated ^ pca k- However, this is not a very likely expla- 
nation. First, as mentioned in § 4.2, D'Elia, Padovani, & 
Landt (2003) have shown that the disk (thermal) compo- 
nent in DXRBS FSRQ is only ~ 15% on average. Second, 
we find no correlation between Lblr or ^ pca k and the ratio 
of disk to total emission as defined in D'Elia, Padovani, & 
Landt (2003). This would be expected if larger BLR lu- 
minosities and/or larger peak frequencies were due to a 
stronger disk component. 

It is more difficult to estimate the total kinetic jet power 
Ljot as it depends on many uncertain astrophysical param- 
eters. We have derived it according to the prescriptions 
of D'Elia, Padovani, & Landt (2003), which we summa- 
rize here briefly. We have used a theoretical relationship 
obtained by Willott et al. (1999) to link jet power to ex- 
tended radio emission at 151 MHz (using c = 20 in cq. 10 
of D'Elia, Padovani, & Landt (2003), which agrees with 
another independent method), estimating the latter from 
the total power at 5 GHz and a correlation between core- 
dominance parameter and radio spectral index. We stress 
that, although the precise values of Lj ot for a given source 
can be somewhat uncertain, we believe that its statistical 
use is warranted. 

Figure 9 plots the kinetic jet power vs. ^ poa k for our 
sources. No correlation is present between the two quanti- 
ties. As Lj e t is an intrinsic power, again this is in contrast 
with the "blazar sequence" of Fossati et al. (1998) and 
Ghiscllini ct al. (1998). Similarly to Fig. 7, however, in 
Fig. 9 there might be indications of an upper envelope in 
the right part of the diagram. For example, sources with 
Ljot > 10 39 W all have ^ poa k ^ 10 15 Hz while sources with 
L r > 10 40 W have f peak < 10 14 Hz. 



Fig. 8. — Broad Line Region luminosity vs. the synchrotron peak 
frequency for the DXRBS FSRQ. See text for details. 
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Fig. 9. — Kinetic jet power vs. the synchrotron peak frequency 
for the DXRBS FSRQ (filled points) and BL Lacs (crosses). See 
text for details. 



tral fit to the full pulse-height analyzer (PHA) spectrum. 
One important limitation of this method is the effect of 
the PSPC background, which is not subtracted off the 
WGACAT hardness ratios and becomes more and more 
important close to the sensitivity limit, especially at large 
off-axis angles. Fiore et al. (1998) have discussed this at 
length and shown that for sources with signal-to-noise ra- 
tio (SNR) > 7 background contamination is not impor- 
tant. However, many of the Fiore et al. (1998) sources 
were targets and therefore at relatively small PSPC off- 
sets. As DXRBS is a serendipitous survey, we have ex- 
cluded all targets and therefore, on average, our sources 
are at larger PSPC offsets, where the PSPC point spread 
function (PSF) degrades significantly, increasing the seri- 
ousness of the background contamination. We have then 
conservatively chosen a higher SNR cut of 10, which re- 
duces our sample to 54 sources, 39 FSRQ and 15 BL Lacs. 
We stress that these effective spectral indices should still 
be regarded as an estimate of the "average" X-ray spectral 
shape and are therefore most suitable for statistical stud- 
ies. Errors on these spectral indices (la) were derived as 
described in Padovani & Giommi (1996) and are typically 
- 0.2. 



6. FURTHER CLUES 

We have shown the existence of significant numbers of 
FSRQ with broad-band spectra similar to those of HBL, 
that is with synchrotron peak frequencies higher than 
those of "classical" FSRQ and reaching the UV band. We 
now concentrate on DXRBS to further investigate the im- 
plications of our findings. 

6.1. X-ray Spectral Slopes 

The study of the X-ray emission of the FSRQ with syn- 
chrotron peak frequencies in the UV/X-ray band is vi- 
tal. The X-ray emission of HFSRQ, in fact, should be 
synchrotron in nature, in contrast to that of most FSRQ, 
where it is dominated by inverse-Compton emission (e.g., 
Sambruna et al. 1996); a dichotomy similar to that exhib- 
ited by LBL and HBL. In BL Lacs, in fact, observational 
evidence points to a different origin for the X-ray emis- 
sion of HBL and LBL (e.g., Padovani & Giommi 1996). 
In HBL, the X-ray continuum appears to be an exten- 
sion of the synchrotron emission seen at lower energies, 
consistent with their steep (a x ~ 1.5) X-ray spectra in 
the ROSAT band. In LBL, the X-ray continuum is more 
likely due to inverse Compton emission, consistent with 
their harder (a x ~ 1) spectra. BeppoSAX observations of 
BL Lacs are confirming this picture (Wolter et al. 1998; 
Padovani et al. 2001). Indeed, Padovani et al. (2002) have 
presented evidence for at least one HFSRQ with relatively 
high 

Vpeak ~ 2x 10 16 Hz and steep (a x ~ 1.5) synchrotron 
X-ray spectrum and two possible "intermediate" sources 
with v pcak « 10 15 Hz. 

We have derived X-ray spectral indices in the 0.4 — 2.0 
keV range for the DXRBS FSRQ using hardness ratios 
following Padovani & Giommi (1996). This method, when 
applied to relatively bright X-ray sources, is relatively ro- 
bust and compares well with the results of a proper spec- 
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Fig. 10. — ROSAT X-ray spectral index vs. the synchrotron peak 
frequency for DXRBS FSRQ (filled points) and BL Lacs (crosses). 
Error bars represent la uncertainties. 

Figure 10 shows the ROSAT X-ray spectral index vs. 
the synchrotron peak frequency for DXRBS FSRQ (filled 
points) and BL Lacs (crosses). We note that a x is rela- 
tively flat (~ 1 — 1.5) and constant for ^ pca k ^ 10 14 Hz. 
For 10 14 < ^pcak ^ 10 15 Hz, a x steepens to reach values 
up to ~ 2.5. Above f p0 ak ~ 10 15 Hz a x flattens again 
to reach values ~ 1 — 1.5. Fig. 10 shows a trend similar 
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to that displayed by the BL Lacs included in Fig. 6 of 
Padovani & Giommi (1996), despite the fact that ~ 70% 
of the sources are FSRQ. We then infer that the interpreta- 
tion put forward for BL Lacs applies to our FSRQ as well. 
Namely, at low f pca k values flat inverse Compton emission 
dominates. For "intermediate" values the steep tail of the 
synchrotron component enters the ROSAT band. Finally, 
when ^p ea k gets even closer to the X-ray band, the X-ray 
spectrum will flatten out again, because the ROSAT band 
is now sampling the top of the synchrotron emission. 
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Fig. 11. — Radio-X-ray spectral index ct rx vs. synchrotron peak 
frequency fpeak f° r FSRQ (filled points) and BL Lacs (crosses) 
for the DXRBS sample. The dotted lines at a TX = 0.725 and 
logi/p Cak = 15.3 denote the two quadrants (top-left and bottom- 
right) occupied by the sources studied by Fossati et al. (1998). See 
text for details. 

We note that BeppoSAX observations of four HFSRQ 
candidates have given mixed results: one source is clearly 
synchrotron dominated, another one is clearly inverse 
Compton dominated, while two others have a flat X-ray 
spectrum with evidence of steepening at low energies, sim- 
ilar to intermediate BL Lacs (Padovani et al. 2002). It 
is, however, worth noting that the synchrotron-dominated 
source is the one with the lowest a rx and right in the mid- 
dle of the HBL box, while the Compton-dominated ob- 
ject has the largest a rx and is at the edge of the box. 
BeppoSAX spectra of "classical" (2 Jy) FSRQ, on the 
other hand, are all consistently flat, with a x ~ 0.7, per- 
fectly explained by inverse Compton emission (Tavecchio 
et al. 2002; Giommi et al. 2002a). Sambruna, Chou, & 
Urry (2000) reported on ASCA observations of four FSRQ 
characterized by steep ROSAT spectra (a x ~ 1-3). The 
sources were all found to have flat hard X-ray spectra, 
with a x ~ 0.8. Sambruna et al. discuss their results in 
terms of relatively high synchrotron peaks and thermal 
emission extending into the X-ray band. Importantly, 
their sources sample a region of parameter space widely 
different from ours. Their effective spectral indices place 



them firmly in the LBL region, unlike ours, so that they 
should not have been expected to show high ^ pca k values 
and steep ASCA spectra. 

As shown in § 4.2, the f pea k distributions are continu- 
ous and peak at w 10 14 Hz for both FSRQ and BL Lacs. 
While on this basis there appears to be no need for divid- 
ing blazars into "H" (high-energy) and "L" (low-energy) 
subclasses, the picture that comes out for FSRQ is the 
same as that already known for BL Lacs (e.g., Padovani 
& Giommi (1996)), that is one of a division which can be 
based on physical grounds with "H" sources as those whose 
X-ray band is dominated by synchrotron emission and "L" 
sources as those in which inverse Compton dominates. 

6.2. Synchrotron Peak Frequency vs. Effective Spectral 

Indices 

It has been suggested in the literature (e.g., Padovani & 
Giommi 1995b, 1996; Fossati et al. 1998) that the blazar 
synchrotron peak frequency can be estimated from the val- 
ues of the effective spectral indices. We address this point 
here by using for the first time a large, homogeneous and 
well-defined sample of blazars. 

Figure 11 plots a rx vs. ^ poa k- As can be seen, the two 
parameters correlate quite well (the correlation is signif- 
icant at the > 99.9% level for the whole sample and the 
FSRQ and at the 98% level for BL Lacs). However, there is 
considerable scatter - for any given value of a rx the scatter 
in z/peak is more than a decade. This is much more than re- 
ported in this relationship by Fossati et al. (1998), whose 
Figure 8 shows a correlation with typically less than a 
decade of scatter, particularly for f poa k > 10 14 Hz. A care- 
ful comparison shows that the major difference between 
the two plots occurs at low a lx . The plot in Fossati et al. 
shows that for ^ poa k ^ 10 15 3 Hz there are no objects with 
a rx < 0.75, while for f poa k > 10 15 3 Hz the exact opposite 
is the case. We have drawn these regions on Figure 1 1 (dot- 
ted horizontal and vertical lines). With DXRBS, however, 
we see a different picture. As Figure 11 shows, DXRBS 
has effectively expanded the accessible region of parameter 
space by making accessible the region 0.65 < a IX < 0.75 
in the bottom left quadrant. The resulting increase in ac- 
cessible parameter space is huge: while for objects with 
^poak ^5 10 13 Hz this is only about 20%, it is more 
than 50% for objects with 10 14 < ^ poa k < 10 15 Hz. A 
look at Figure 11 shows that while the correlation is still 
highly significant, the plot looks very different at the high 
a rx end than it does at low a rx . At high a rx we see a 
clear envelope, while it is obvious that we do not see an 
envelope at low a Ix . In other words a DXRBS blazar with 
a IX ~ 0.7 is just as likely to have ^ poa k ~ 10 12 Hz as it is 
to have f pea k ~ 10 16 Hz. We note that the sources with 
relatively low a rx and f pea k values are the DXRBS sources 
in the upper-left part of Fig. 1 which have low a IX but are 
outside of the HBL box. 

What is the reason behind this difference? In the 
first place the models arc different: Fossati et al. (1998) 
used a simple third-degree polynomial fit, while the model 
we use is based upon a calculation of synchrotron and 
inverse-Compton spectra for reasonable physical param- 
eters. More importantly, however, all the samples used 
by Fossati et al. - the 1 Jy, 2 Jy and Slew surveys - are 
classic, high flux limit surveys in the radio and X-rays re- 
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spectively. As shown in § 3, this type of survey has the 
effect of imposing what amounts to severe selection biases. 
This can in fact be seen in Figure 8 of Fossati et al. (1998), 
where both their ;/ pea k — a rx arid f pC ak — a ro plots are seen 
to be discontinuous, hardly surprising considering that all 
but one of the objects with ^ pca k > 10 153 Hz in the Fos- 
sati et al. study is found in the Einstein Slew Survey but 
not the 1 Jy (only 2/34 1 Jy BL Lacs fall in the lower 
right hand quadrant of our Figure 11 as denoted by the 
two dotted lines, compared to 50/60 Slew BL Lacs), while 
the vast majority of the objects with f pea k < 10 15 3 Hz 
come from the 1 Jansky (33/34 1 Jy BL Lacs are in this 
range compared to 10/60 Slew BL Lacs). 
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Fig. 12. — Radio-optical spectral index a IO vs. synchrotron peak 
frequency ^ pea it for FSRQ (filled points) and BL Lacs (crosses) for 
the DXRBS sample. The dotted line at a ro = 0.825 denotes the 
limit of the sources studied by Fossati et al. (1998). Sec text for 
details. 

A tighter correlation appears to be present between a ro 
and ^pcak, shown in Fig. 12. The correlation is significant 
at the > 99.99% level (~ 97% for BL Lacs only), with 
c^ro oc Vpegk 2 - This correlation is expected to break down 
for t'peak 10 16 Hz or a ro < 0.4 (Padovani & Giommi 
1995b; Fossati et al. 1998) but most of our sources are 
outside of this range, which explains why the correlation 
looks almost linear. Apart from the lack of objects with 
a m ^ 0.8 and the paucity of objects with intermediate 
values of f pca k Fig. 8 of Fossati et al. (1998) docs not 
look too different from Fig. 12 in the overlapping range of 

^peak • 

We notice that Fig. 11 shows that a definition based 
solely on X-ray-to-radio flux ratio or a rx is not optimal, as 
we have found blazars with low a rx and low f poa k values. 
The position of the sources on the a OXl a ro plane, which 
means using two (instead of one) effective spectral indices, 
on the other hand, appears to be more sensitive to the syn- 
chrotron peak frequency, with a difference of a factor ~ 35 
in the mean ^ pca k values for blazars in and out of the HBL 



box. This is true especially for FSRQ. Fig. 1 shows, in 
fact, that most of the blazars with a rx < 0.78 but outside 
the HBL box are strong-lined. Tab. 1 shows also that, 
taking into account the sky coverage, while ~ 54% of BL 
Lacs with a rx < 0.78 are in the HBL box, this is true only 
for - 36% of the FSRQ. 

6.3. Other Properties of HBL/LBL and HFSRQ/LFSRQ 

We have studied the properties of blazars in and out of 
the HBL box not addressed in § 4, namely radio, optical, 
and X-ray powers and redshift. As before, distributions 
are "corrected" by using the sky coverage. HFSRQ have 
slightly smaller radio powers and slightly larger X-ray pow- 
ers than LFSRQ, by factors ~ 2, but not significantly so 
(P ~ 90%). HFSRQ, on the other hand, have significantly 
(P > 99.99%) larger (factor <~ 7) optical powers, due to 
the fact that for these sources (t'peak) ~ 7 x 10 14 Hz. Both 
classes have similar mean redshifts, ~ 1.6 and ~ 1.7 re- 
spectively, with redshift distributions which arc different 
at the 96% level. The same comparison in the case of BL 
Lacs is hampered by the relatively small number statistics 
(18 LBL and 4 HBL with redshift). With this caveat, we 
find that HBL have slightly larger radio and optical pow- 
ers than LBL, by factors ^2 — 5, but not significantly so 
(P ~ 50% and ~ 90% respectively), while they have sig- 
nificantly (P ~ 98.7%) larger (factor <~ 35) X-ray powers. 
The redshift distributions for the two classes are similar 
with (z) ~ 0.4. 

If indeed all blazars that lie within the HBL "box" have 
radio-to-X-ray continua produced by synchrotron emis- 
sion, one might expect a significant population of HFSRQ 
to be 100 GeV - TeV emitters. An example of such an 
object might be RGB J1629+4008 (Padovani ct al. 2002), 
which has a synchrotron peak ~ 2 x 10 16 Hz and a pre- 
dicted vF v at 10 25 Hz of nearly 10~ 10 erg cm~ 2 s -1 . A 
similar point was also made by Perlman (2000), albeit 
without the aid of BeppoSAX spectral data, who also 
listed five FSRQ with a rx < 0.78 and z < 0.1 derived 
from the RGB and Einstein Slew Survey. Thus the dis- 
covery of a new, X-ray loud population of FSRQ predicts 
the existence of a new class of 100 GeV - TeV emitting 
sources which could be particularly helpful for probing the 
diffuse infrared background at higher redshifts. 

7. SUMMARY AND CONCLUSIONS 

We have used the results of two recent surveys, DXRBS 
and RGB, to study the spectral energy distribution of 
about 500 blazars. Never before had this been done with a 
sample even remotely close to ours in terms of size, depth, 
and well-defined selection criteria. DXRBS, in particular, 
is ~ 95% complete and reaches fluxes ~ 20 times lower 
than previously available blazar radio surveys. We have 
first derived the effective spectral indices a ox , a ro , and a rx . 
We have then studied the a rx distributions for DXRBS, 
RGB, and even the EMSS samples for both FSRQ and 
BL Lacs. The serendipitous nature of DXRBS and the 
EMSS has been taken into account by "correcting" these 
and other distributions using the appropriate sky coverage. 
The synchrotron peak frequencies for DXRBS blazars have 
also been derived by using multi-frequency information 
and an homogeneous synchrotron - inverse self-Compton 
model. Broad line region and jet powers were also esti- 
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mated. One of the main aims of this work was to look for 
the strong- lined counterparts of high-energy peaked (HBL) 
BL Lacs, the HFSRQ, that is FSRQ with high-energy syn- 
chrotron peaks. We have found them. Our main results 
can be summarized as follows: 

1. About 10% of DXRBS FSRQ have effective spec- 
tral indices typical of HBL (to be compared with 
15% for BL Lacs) and can therefore be called HF- 
SRQ. The fractions of HFSRQ and HBL increase 
to - 30% and - 60% for RGB and to - 80% and 
100% for the EMSS, respectively. Although HF- 
SRQ have X-ray-to-radio flux ratios larger than pre- 
viously known FSRQ, in none of the samples they 
manage to reach values as high as those of HBL. 

2. The synchrotron peak frequency distribution of 
DXRBS FSRQ and BL Lacs is continuous and 
peaks at « 10 14 Hz, with the former sources hav- 
ing an average ^ pca k ~ 2.5 times smaller than the 
latter. We have verified that blazars with effective 
spectral indices typical of HBL indeed have larger 
t'poak values (by a factor ~ 35) than other blazars. 
About 60% and > 8% of DXRBS BL Lacs have 
^peak > 10 14 and 10 15 Hz respectively, to be com- 
pared with - 30% and - 5% for FSRQ. 

3. These results, together with the dependence we find 
of the X-ray spectral index, estimated from the 
hardness ratios, on f poa k, confirm the existence of 
strong-lined counterparts of high-energy peaked BL 
Lacs. As is the case for HBL, we would expect a 
significant fraction of these sources to emit at 100 
GeV - TeV energies. 

4. We find no anti-correlation between synchrotron 
peak frequency and radio, broad line region, and 
jet powers, contrary to the predictions of the so- 
called "blazar sequence" scenario, which calls for 
an inverse dependence of ^ pca k on intrinsic power 
due to the effects of the more severe electron cool- 
ing in more powerful sources. On the other hand, 
available data from DXRBS and other surveys sug- 
gest that high-^p 0a k-high-power blazars have not 
been found yet, and that HFSRQ do not reach the 
extreme synchrotron peak frequencies of BL Lacs. 
This indicates that after all there might be an in- 
trinsic, physical limit to the synchrotron peak fre- 



quencies and therefore electron energies which can 
be reached by powerful blazars. 

The discovery of HFSRQ and the study of their proper- 
ties have important implications for our understanding of 
jet formation and physics. Since ^ pC ak °c 7p ea k^' wnere 
7peak is the Lorentz factor of the electrons emitting most of 
the radiation, we have shown that powerful jets with large 
magnetic fields and electron Lorentz factors can indeed ex- 
ist - regardless of whether or not they have strong emission 
lines - albeit up to a point. This provides an important 
challenge for existing models that advocate that the spec- 
tral energy distribution of relativistic jets is strongly af- 
fected by the external radiation field. We have also shown 
that selection effects are very strong and that, in particu- 
lar, the HBL/HFSRQ fraction is sample-dependent. Can 
we separate selection effects from physics? We think we 
can. By using DXRBS we have sampled a region of pa- 
rameter space which should be largely unbiased in terms 
of fpeak (unlike that covered by the EMSS, for example). 
We have found that HBL/HFSRQ make up a minority of 
the blazar population, ~ 10 — 15%. We then believe that 
the available evidence suggests that nature preferentially 
makes jets which peak at IR/optical energies. We will 
address this issue in detail in a future paper. 

It is also clear that, although a consistent picture comes 
out of our results, we need more information on these 
sources. We are planning dedicated X-ray observations 
of our newly discovered HFSRQ to further constrain their 
X-ray emission processes. 
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